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Novel Energy Scale in the Interacting 2D Electron System Evidenced from Transport 
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By analyzing the in-plane field magnetoconductivity, zero field transport, and thermodynamic spin 
magnetization in 2D correlated electron system in high mobility Si-MOS samples, we have revealed 
a novel high energy scale T*, beyond the Fermi energy. In magnetoconductivity, we found a sharp 
onset of the novel regime Sa(B,T) oc ( B/T ) 2 above a density dependent temperature Tkmk(n), 
the high-energy behavior that “mimics” the low-temperature diffusive interaction regime. The 
zero field resistivity temperature dependence exhibits an inflection point Tinfl. In thermodynamic 
magnetization, the weak field spin susceptibility per electron, d\/dn changes sign at TdM/dn■ All 
three notable temperatures, Tkink, Tin a, and TdM/dn , behave critically oc (n — n c ), are close to each 
other, and are intrinsic to high mobility samples only; we therefore associate them with a novel 
energy scale T* caused by interactions in the 2DE system. 

PACS numbers: 71.30.+h, 73.40.Qv, 71.27.+a 


Two-dimensional (2D) interacting low density carrier 
systems in the past two decades attracted consider¬ 
able interest [l, 3], demonstrating fascinating electron- 
electron interaction effects^ such as metallic temperature 

metal-insulator transition 


dependence of resistivity c _ 

(MIT) @,0, strong positive magnetoresistance (MR) in 
parallel field , strong renormalization of the effective 
mass and spin susceptibility 0,0 , etc. 

Far away of the critical MIT density n c , in the well 
“metallic regime”, these effects are explained within 
framework of the Fermi liquid theory - either in terms 
of interaction quantum corrections (IC) pm , or tem- 
rerature dependent screening of the disorder potential 


peraturc 

mm. 


Both theoretical approaches so far are used to 
treat the experimental data on transport, and the former 
one - also to determine the Fermi liquid coupling con¬ 
stants from fitting the transport and magnetotransport 
data to the IC theory. In the vicinity of the critical re¬ 
gion, conduction is treated within renormalization group 
17], or Wigner-Mott approach 18, [l9|- This regime is 


however out of scope of the current paper. 

On the other side, numerous theories predict break¬ 
down of the uniform paramagnetic 2D Fermi liquid state 
as interaction strength increases 2Cj-24|. However, it re¬ 
mained almost unexplored how the potential instabilities 
may reveal themselves in thermodynamics and transport. 

In this paper we report results of the transport, mag¬ 
netotransport and magnetization measurements with 2D 
correlated electron system, which reveal the existence of 
a novel characteristic energy scale T*, that is smaller 
than the Fermi temperature Tp, but much bigger than 
1/r (we set throughout the paper h,kB,PB — !)• Ob¬ 
viously, no such large energy scale may exist in the 
pure Fermi liquid. T* reveals itself (i) in the weak in¬ 
plane field magnetotransport, (ii) in zero field trans¬ 
port, and (iii) in the spin magnetization per electron. 
In magnetoconductivity, we found a sharp onset of the 


novel regime Sa(B,T ) oc (B/T) 2 above a density de¬ 
pendent Tkink (n), the high-energy behavior that “mim¬ 
ics” the low-temperature diffusive interaction regime [jjJ . 
Tkink (n) correlates well with inflection point Tinfl (n) in 
the zero field resistivity temperature dependence. Fi¬ 
nally, the two remarkable temperatures correlate with 
the temperature TdM/dn for which the spin susceptibil¬ 
ity per electron dx/dn (and dM/dn) changes sign. All 
three notable temperatures, Tkink, Ti n fl, and T dM / dn , be¬ 
have critically oc (n — n c ), are intrinsic to high mobility 
samples only, and are close to each other; we therefore 
associate them with a novel energy scale T* caused by 
interactions in the 2DE system. 

Experimental. The ac-measurements (5 to 17 Hz) of 
resistivity were performed using the four-probe lock-in 
technique in magnetic fields up to ±7T. The range of 
temperatures, 0.4 — 20 K, was chosen so as to ensure the 
absence of the shunting conduction of bulk Si at the high¬ 
est temperatures, and, on the low-temperature side, to 
exceed the valley splitting and intervalley scattering rate 
[25l ] . The studied high mobility (lOO)Si-MOS samples 
had w 190 nm gate oxide thickness, and were lithograph¬ 
ically defined as rectangular Hall bars, 0.8 x 5 mm 2 [201 ]. 
The magnetoconductivity measurements are performed 
similar to Ref. [27| , but in the much wider domain of 
densities and temperatures, from far above the MIT crit¬ 
ical density and in the well-conducting regime kpl 3> 1 
down to the critical regime kpl ~ 1. 

By rotating the sample with a step motor, we aligned 
magnetic field in the 2D plane to within 1' accuracy, 
using the weak localization magnetoresistance as a sen¬ 
sor of the perpendicular field component. Carrier den¬ 
sity n was varied by the gate voltage V g in the range 
(0.9 — 10) x 10 11 cm” 2 . The linear n(V g ) dependence was 
determined from quantum oscillation period measured 
in the perpendicular held orientation during the same 
cooldown. 
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In-plane field magnetoconductivity. The inset in 
Fig. Q] shows that the magnetoresistivity varies in weak 
fields as B 2 , with a high accuracy. From the p{B) 
data we determined the magnetoconductivity prefactor 
a a = —\d 2 a I 8B 2 \b=q = (l/2p 2 )d 2 p/dB 2 \ B =o which is 
analyzed below. 



T (K) 

FIG. 1: (Color online) Temperature dependence of resis¬ 
tivity at zero field in the vicinity of n c (« 0.85) for sam¬ 
ple Si-2. The densities are in 10 11 cm -2 . Crosses mark 
the p(T) maxima, dots - the inflection points. The inset 
shows p versus B 2 for sample Si-4 at a fixed density and 
at five temperatures: 0.4, 0.9, 1.3, 1.8, and 2.5 K (from 
top to bottom). Vertical ticks mark B = T field limit. 


In the weak field limit, B < T, variations of the con¬ 
ductivity at a fixed temperature are low, < 5%, (see 
insert to Fig. 1). This smallness favors comparison of 
the data with theory of interaction corrections which 
makes firm predictions specifically for magnetoconduc¬ 
tivity (MC) EH- The magnetoconductivity prefactor 
a a (T,n) is plotted in Fig. [2] versus temperature. In the 
wide density range, (2 — 10) x 10 11 cm -2 in the well- 
conducting regime, the estimated diffusive/ballistic bor¬ 
der Em , T db = (1 + Fq)/2ttt « 0.2 K, is below the ac¬ 
cessible temperatures range of our measurements and we 
anticipate to observe the ballistic regime of interactions. 

Surprisingly, at temperatures much higher than Tdb, 
the prefactor a a (T) develops roughly oc T~ 2 (contrary 
to the predicted ballistic-type dependence a a oc T _1 ) 
EH- At somewhat lower temperatures (but still higher 
than Tdb), the a a (T) dependence softens to oc T _1 . The 
crossover in Fig. [2] occurs rather sharply, as a kink on 
the double-log scale. The kink and the overall type of 
behavior is observed in the wide range of densities and 
for all studied high mobility samples. Figure 3 shows 
that Tkink(?T-) develops critically versus electron density, 
oc (n — n c ) where n c within experimental uncertainty co¬ 
incides with MIT critical density. The sharp crossover at 
high temperatures to the novel regime of MC, which is in 


contrast with the theory predictions, is one of the main 
results of our study. (In [281 ] we compare the magnetore¬ 
sistivity and the magnetoresistivity prefactors). 



FIG. 2: (Color online) Temperature dependence of the 
magnetoconductivity prefactor a a for sample Si-2, for 
eight electron densities (increasing from bottom to top), 
in units of 10 n cm~ 2 . For clarity, the curves are mag¬ 
nified by the factors shown next to each curve. Arrows 
mark the kink positions, the dashed curves show Tp(n) 
and Tp(n), dash-dotted curve shows Ti nco h(n). 

Features in thermodynamics. In case the kink in 
magnetoconductivity indeed signals a novel energy scale, 
it must show up in temperature dependencies of other 
physical quantities measured in the high temperature 
range and in weak or zero magnetic fields. Other avail¬ 
able data which fits these requirements are as follows: (i) 
spin magnetization per electron dM/dn, and (ii) zero- 
field transport p(T). 

The spin magnetization data [29 ] show a pronounced 
sign change of dx/dn = d 2 M/dBdn at a density depen¬ 
dent temperature T d M/dn(n)- Physically, the sign change 
means that for temperatures lower than T dM / dn (n), the 
minority phase (large spin collective “spin droplets”) 
melt as density increases. In other words, extra electrons 
added to the system join the Fermi sea, improve screen¬ 
ing and favor “spin droplets” melting. For temperatures 
above T d M/dn{n ), the number of “spin droplets” grows as 
density increases; here the extra electrons added to the 
2D system prefer joining the “spin droplets” (see also 
[H]). The T d M/ d r,(n) dependence copied from Figs. 1 
and 2 of Ref. [29| is depicted in the insert to Fig. [3] 
One can see that T d M/dn{n) behaves critically and van¬ 
ishes to zero at n c ; remarkably, within the measurements 
uncertainty, it is consistent with Tki n k(n) deduced from 
magnetotransport. 

Zero field transport. Figure [Q shows typical tem¬ 
perature dependencies of resistivity for high mobility Si- 
MOS samples. Each curve has two remarkable points: 
the p(T) maximum, T max , and inflection, Ti n fl .3Cj. 
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Whereas T max is an order of the renormalized Fermi en¬ 
ergy, the inflection point happens at much lower temper¬ 
atures, in the degenerate regime. Importantly, the inflec¬ 
tion temperature appears to be close to the kink tempera¬ 
ture (see Figs. mm- Therefore, the proximity of the three 
notable temperatures which are inherent to high mobility 
samples solely, T kink ss T infl « T dM/dn strongly suggests 
the existence of a new energy scale T* in the correlated 
2D system. T* is much less than the bare Fermi temper¬ 
ature Tp [3lj |. and the renormalized Tp = Tp(rrib/m*) 
[0 ; in contrast to Tp (which is oc n), T*(ri) develops as 
(n — n c ). On the other hand, it is much higher than the 
“incoherence” temperature at which the phase coherence 
is lost (defined as t v (T) = r 0 ], confirming that the 
kink, inflection and dx/dn sign change are irrelevant to 
the single-particle coherent effects. 



FIG. 3: Empirical phase diagram of the 2DE system. 
Dashed areas are: (I) - the ballistic interaction regime, 
(II) - novel MC regime. Hatched area (III) is the non¬ 
degenerate regime, blank area localized phase. Full 
dots - the kink temperature T k j nk , open dots - inflec¬ 
tion point Ti n fl. Sample Si-2. Dashed curves show the 
calculated bare (Tp) and the renormalized (Tp) Fermi 
temperatures. Insert blows up the low density region; 
dashed line is T dM / dn 0 ], 

Phenomenological model for transport and 
magnetotransport. In the absence of an adequate mi¬ 
croscopic theory, we attempt to elucidate the origin of 
the T* energy scale and of the novel magnetoconductance 
behavior. We suggest below a phenomenological model 
that links ‘“high temperature” transport and magneto¬ 
transport behavior in a unified picture. 

The features of our interest, T k j nk and Ti n fl represent 
“high-energy” physics. Moreover, the p(T) (and a(T)) 
variations of the experimental data (Fig. [TJ are so large, 
that the first order in T corrections, of cause, cannot 
describe them. Our analysis of other known theoretical 
models for homogeneous 2D Fermi liquid 0] (see also 
0 ) reveals that neither of them describes adequately 


the inflection on the p(T) data and of cause does not 
include an associated energy scale. For this reason, we 
turn attention to the two-phase state. 

There is a large body of theoretical suggestions for 
spontaneous formation of the two phase state 00,0 
24 . 3l| due to instabilities in the charge or spin channel, 
and many experimental indications obtained with meso¬ 
scopic systems or local probes 39, [37]]. Finally, the spin 
magnetization measurements |29( revealed the existence 
of the two-phase state of the macroscopic 2DE system, 
where the minority phase - large spin droplets - coexist 
with the majority Fermi liquid state. Dealing with the 
two-phase state, the two channel scattering or additive 
resistivity approach seems quite adequate to the problem. 

The typical functional form of the p(T) (Fig. 1) also 
prompts dual channel scattering. The simplest functional 
dependence that correctly describes the inflection in p(T) 
is provided by the phenomenological form Si: 


p(T) = po + pi exp - 


A (n) 
~T~ 


A (n) = a(n - n c (B)), 

(1) 

where p\(n,B) is a slowly decaying function of n, and 
po(n,T) includes Drude resistivity and quantum correc¬ 
tions, both from the single-particle interference and in¬ 
teraction. Although this model was suggested on a dif¬ 
ferent ground, it fits well resistivity data in the vicinity 
of MIT for various material systems (see 01) ■ This sim¬ 
ple additive p(T) form satisfies general requirements for 
the transport behavior in the vicinity of a critical point 
130 . 0, and also explains the apparent success of the 
earlier attempts of one-parameter scaling (namely of the 
p(T) steep rise, mirror reflection, etc.) )3|, |5|. 

Obviously, in this model Ti„fl = A/2. To take magnetic 
field into account, and following results of Refs. 0, 0] 
we include to (A/T) all the lowest order in B/T (and 
even-in-H) terms, as follows: 


A(T, 5, n)/T = A 0 (n)/T-P(n)B 2 /T-an)B 2 /T 2 , (2) 
with Ao = a[n — n c (0)]. 

Eqs. ([T|) and (|2|) link the magnetoconductance with the 
zero-field p(T) temperature dependence. With these, the 
p(T , B ) dependence is as follows: 


p(B,T ) = [& D 5<j ■ exp ( T/Tb)) 

n - n e (0) _ 
T 


+ pi exp —a 


-1 

B 2 B 2 

@ T ^ J^2 


.(3) 


The term in the square brackets includes the Drude 
conductivity and interaction quantum corrections [11 . 
IT2I ]. The latter, Sa(T) = j(B 2 /T) + 77T, was calculated 
using experimentally determined Ftf (n) values 0, 0, 
and o d found from a standard procedure 0 - In order 
to cut-off the corrections above a certain border temper¬ 
ature (0 and, thus, to disentangle the exponential- and 
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linear-in-T contributions, the calculated interaction cor¬ 
rection are cut-off with an exponential crossover function 
above Xb which for simplicity we set equal to A(n)/2. 

From Eq. (ED, the prefactor a a = —(l/2)d 2 a/dB 2 is 
calculated straightforward and in Fig. [4] is compared with 
experimental data. In the p(T) fitting [Figs. 4(a,c,e,g)], 
basically, there is only one adjustable parameter, pi(n), 
for each density. Indeed, n c (0) is determined from the 
conventional scaling analysis at B = 0 [30], and the slope, 
a = 2dT- ln fi(n)/dn may be determined from Fig. [3] How¬ 
ever, in order to test the assumed linear A (n) relation¬ 
ship, Eq. (1), we treated a(n) as an adjustable parameter. 
On the next step, in the a a (T) fitting [Figs. 4 (b,d,f,h)] 
we fixed the parameters determined from the p(T) fit, 
and varied /3(n) and £(n). 



FIG. 4: Fitting p(T, B = 0) dependencies (left panels) 
and a a (T) (right panels) with the same set of the fit¬ 
ting parameters. Sample Si-2; carrier densities (from top 
to bottom) are n = 1.5; 2.0; 2.5, and 3.25 x 10 11 cm -2 . 
Fitting parameters are presented in the Table. 


One can see that both p(T) and a a (T) are well fitted; 
the model captures correctly the major data features, the 
steep p(T) rise (including the inflection), and the a„(T) 
kink. Within this model, the kink signifies a transition 
from the low-temperature magnetoconductance regime of 
ballistic interaction (where the exponential term may be 
neglected and where the quantum corrections dominate) 
to the high temperature regime governed by the steep 
exponential p{T) rise. The parameters of the fit (Figs. [4]) 
are summarized in the Table of the Supplementary mate¬ 
rials. The factor /3 is an order of magnitude smaller than 
£, therefore, the corresponding term in Eq. (3) becomes 
important only at high temperatures. The slope, a, is 
almost constant, confirming our assumption [(Eq. (1)]. 


Impact on the MC interpretation. For high den¬ 
sities n n c , the temperature range above Tkmk is un¬ 
ambiguously beyond the diffusive regime of interactions 
and, hence, the B 2 /T 2 dependence is the novel high- 
temperature MC regime of the non-diffusive type. Below 
Xkink the temperature is still higher than Tdb and the 
regime a a oc T^ 1 (see Fig. [2J therefore is reminiscent of 
the standard ballistic interaction regime [12]. This con¬ 
clusion is confirmed by Figs. 4b,d,f,h where the standard 
interaction corrections (incorporated in Eq. 3) with ex¬ 
perimentally determined interaction parameters provide 
quite a sucsessful fit below Tkink- 

The two-stage “high-temperature” dependence of the 
magnetoconductivity prefactor is universal and persists 
from well conducting regime far above the critical den¬ 
sity, down to the very critical density. As a result, in 
the vicinity of MIT, the “high-temperature” MC regime 
da oc —(B/T) 2 , mimics the behavior anticipated for the 
diffusive regime of quantum interaction [12l.l46j|. This fact 
therefore casts serious doubt on the RG treatment of the 
experimental p(T,B ||) data in the vicinity of MIT, and 
particularly, on the phase diagram of 2D interacting and 
disordered system deduced from fitting the experimental 
data with the RG-theory (47], Q . 

Conclusion. To conclude, we have found a novel high 
energy scale T* to exist in the correlated 2D electron sys¬ 
tem, beyond the Fermi energy. It reveals itself in trans¬ 
port, in-plane field magnetotransport and thermodynam¬ 
ics. All three notable temperatures behave critically, 
oc (n — n c ), and are rather close to each other. For tem¬ 
peratures above the density dependent T*, the in-plane 
field MC crosses over from the conventional ballistic- 
type — (B 2 /T) dependence to the novel —( B 2 /T 2 ) de¬ 
pendence. We suggested phenomenological description 
of the transport and magnetotransport data, based on 
the two phase state (two resistivity channels). 

Our present studies are performed with high mobil¬ 
ity Si-MOS samples which show well pronounced MIT. 
Therefore, the results obtained send a warning to inter¬ 
pretation of the in-plane field MC in the vicinity of MIT, 
particularly within the framework of the RG theory |4f|. 
Secondly, our results explain why the Fermi-liquid pa¬ 
rameters extracted from fitting the measured MC scatter 
significantly in various experiments: indeed, by fitting 
the data in the nominally ballistic regime, one would 
observe a a (and deduce Ftf values) strongly dependent 
on the particular temperature range, above or below the 
kink. 

Clearly, there is a need in a microscopic theory that 
must explain on the same footing all three critical behav¬ 
iors: in the zero field resistivity, in magnetoconductivity, 
and in spin susceptibility per electron. A possible origin 
of the T* scale may be such a structure of collective en¬ 
ergy levels for individual droplets of the minority phase, 
which in analogy with quantum dots may cause features 
simultaneously in thermodynamics and in transport of 
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itinerant electrons. 
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Comments on calculation of a CT from the measured p(B) dependence 


In weak in-plane fields: 


<(B 2 


a = do — a a B 2 


p — Po + o, p B 2 + o(B~), 


(1) 

( 2 ) 


where B is considered to be small as compared with either T, or ( T 2 t) ( h , ks and ps are set to unity throughout 
the paper, for shortness), and by definition 


= -i<9 2 er/&B 2 


B=0 


2 p = ^ 2 p/dB‘ 


B=0 


Consider the relation between a a and the experimentally measured p(B). In purely parallel magnetic field <r = 1 /p. 
Taking the second derivative from both sides of equation m and recalling that (dp/dB )| B=Q = 0, we obtain: 


1 d 2 p 1 
2p 2 dB 2 ~ 


dp_ 

dB 


J B =0 


1 d 2 p 


(3) 


The latter proves Eq. (1) used in our paper to calculate a a from the measured p(B) dependence. 

In the experimental data, purely parabolic p(B) oc B 2 dependence was found to extend with high accuracy far 
above the range of low fields (B < T) which we used for calculating a a . For example, the higher order-in- (B/T) terms 
in Eqs. (ID were less than 0.1% (relative to the S 2 -term) even at B/T = 6.5, and have been neglected therefore for 
B « T. 


Magnetoconductivity and magnetoresistivity 

Variations of the conductivity with weak in-plane field at a fixed temperature are low, < 5%, in the selected range of 
fields B < T (see insert to Fig. 1 of the main text). This smallness enables to apply the theory of interaction corrections 
(IC) which makes firm predictions for the functional temperature dependence of the weak field magnetoconductivity 
(MC) and suggests a clear physical picture behind it [l]. Therefore, in the main text of our paper we analyze namely 
the MC prefactor, a G = —(l/2)d 2 a / dB 2 \ B —0 versus T. In Fig. la here and in the main text, one can see that a a 
develops oc T -1 in the low-temperature part of the ballistic regime ((1 + F//)/2ttt « 0.2K < T), in a qualitative 
agreement with the IC predictions. However, for higher temperatures, still in the ballistic and, presumably, degenerate 
regime T < Tp , there is a sharp kink and the onset of the novel unforseen dependence, a a (T ) oc T~ 2 . 

This effect is observed for several high and moderate mobility Si-MOS samples (see Ref. [26] in the main text), and is 
missing in more disordered (a factor of 10 lower mobility) samples, where a a {T) dependence is qualitatively consistent 
with theory of interaction corrections. Importantly, by contrast to the Fermi energy which in 2D is proportional to 
the carrier density, the novel energy scale develops critically, oc (n — n c ), with slightly sample mobility dependent n c 
values. It vanishes to zero as density approaches the critical value (see Fig. 3 of the main text); these two facts confirm 
the relevance of the electron-electron interaction effects. Another indication of the crucial importance of the electron- 
electron interactions is the fact that the kink in a CT (T) at Tkink and the novel regime of MC at T > Uink are intrinsic 
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only to high mobility samples, where the strongly correlated regime is accessed upon lowering density. For samples 
Si-40 and Si-46 with a factor of 10 - 30 lower mobilities, in the same range of temperatures, the magnetoconductance 
develops in accord with IC theory predictions for the diffusive regime, with no kink. 

The magnetoconductivity in the in-plane field is inequivalent to magnetoresistivity (MR), because variations of 
the conductivity with temperature at zero field are large (a factor of 4 - 7) in high mobility Si-MOSFETs. As a 
result, the a p {T) = (l/2)d 2 p/dT 2 \B=o temperature dependence is different from that for a„{T) (see Fig. lb): it is 
less transparent, being affected by both, the onset of the novel regime in MC and by the strong p(T) (and a(T)) 
dependencies. For higher densities, where the p(T ) variations are relatively weak (see the lowest 4 curves in Fig. lb), 
a p (T) exhibits a maximum that coincides with the kink in a a (T) for densities n = 10,5.25,3.25,2.5 x 10 1:L cm~ 2 . 
For lower densities the maxima in a p (T) get smeared which hampers their quantifying. The simplicity of the a a {T) 
dependence (in comparison with a p (T)) clearly points at the primary role of the magnetoconductivity rather than 
magnetoresistivity for the studied system. 

The kink temperature Tkink in a a lies far away from the bare and renormalized Fermi energy and from (1+Fq)/2ttt « 
0.2 K value Q which are the only known energy scales in the Fermi liquid. We interpret Tkink as a manifestation of 
an additional energy scale, beside the Fermi energy. Obviously, no such energy scale may exist in the pure 2D Fermi 
liquids, and vice versa, its existence indicates a non-Fermi liquid state. 

In Fig. la, one can also see that the magnetoconductivity prefactor exhibits another twist upward, clearly noticeable 
for the four lowest curves (lowest densities). However, this feature occurs at approximately renormalized Tp and is 
likely to signify a transition to non-degenerate regime, which is beyond the scope of our paper. 




FIG. 1: (Color online) Temperature dependencies of (a) a a (T) and (b) a p . Arrows (left panel) and dots (right panel) mark 
Tkink - kink temperature in a a . 


Other available data: spin magnetization 

In order to test whether the kink temperature in magnetoconductivity has a more general significance, we inspected 
temperature dependencies of other physical quantities in the high temperature range and in weak or zero magnetic 
fields. Other available data which fit these requirements are as follows: 

(i) spin magnetization per electron dM/dn Q, and 

(ii) zero-field transport p{T). 

The spin magnetization-per-electron dM/dn data [5j], in general, are interpreted as a clear evidence for the formation 
of a two-phase state, in which the Fermi liquid phase coexists with large spin collective “spin droplets” (the latter 
being presumably collective localized states). These dM/dn data, in particular, show a pronounced sign change 
of dx/dn = d 2 M/dBdn at a density dependent temperature T dM / dn (n). Physically, the positive sign of the spin 
susceptibility per electron, dx/dn , means that for temperatures higher than T d M/dn(i i), the number of “spin droplets” 
grows as density increases. Here the extra electrons added to the 2D system prefer joining the “spin droplets”, increase 
its magnetization M and cause positive sign of dx/dn. 
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FIG. 2: Empirical phase diagram of the 2DE system. Dashed areas are: (I) - ballistic interaction regime, (II) - novel MC 
regime. Hatched area (III) is the non-degenerate regime, blank area - localized phase. Full dots - the kink temperature Jkink, 
open dots - inflection point Ti n fl. Sample Si-2. Dashed curves show the calculated bare (Tf) and the renormalized ( Tp ) Fermi 
temperatures. Insert blows up the low density region; dashed line is T dM /dn ~ sign change of dM/dn [j|. 


For temperatures lower than T dM / dn (n), the number of “spin droplets” (minority phase) decreases as density 
increases, i.e. the “spin droplets” melt as density increases. In other words, extra electrons added to the system join 
the Fermi sea, improve screening and favor spin droplets “melting” that results in the negative sign of dx/dn. The 
T d M/dni n ) dependence copied from Figs. 1 and 2 of Ref. _3j is depicted in the insert to Fig. [2] One can see that 
T dM /dn{ n ) also behaves critically and vanishes to zero at n c ; remarkably, within the measurements uncertainty, it is 
consistent with Tkink(u) deduced from magnetotransport. 


Other available data: resistivity and conductivity in zero Held 

Searching for manifestation of the novel energy scale in zero-field transport we analyze the p(T) and u{T) depen¬ 
dencies at zero field. The variations of these quantities in the relevant temperature range are large (up to a factor of 
7), making the IC theory inapplicable in this “high temperature” regime. One can see from Fig. 3a below that the 
p(T ) temperature dependence is monotonic up to the limits of degeneracy, T = Tp, and follows one and the same 
additive resistivity functional form (Eq. (1) of the main text) over a wide density range: 

p(T) = po +piexp(-A(n)/T), 

A (n) = a(n - n c (B)), (4) 

where p\(n, B) is a slowly decaying function of n, and po(n,T) includes Drude resistivity and quantum corrections, 
both from the single-particle interference and interaction. 

By inspecting available in literature data by several research groups for p{T) in the low density correlated regime, 
one finds that the above empirical additive resistivity form fits well the p(T ) dependence for a number of material 
systems [l-Q. The additive resistivity form (cf. Matthiessen’s rule) naturally corresponds to the two-phase state of 
the low density 2D electronic system. The two-phase state is experimentally revealed in macroscopic magnetization 
measurements Q, local probe experiments [()}, and also widely discussed in theory fTtHD!]. This empirical additive 
p(T) form satisfies general requirements for the transport behavior in the vicinity of a critical point mm , and 
explains the apparent success of the earlier attempts of one-parameter scaling (namely of the p(T) steep rise and the 
mirror reflection symmetry between p(T) and u(T) on the metallic and insulating side of the MIT) Il8l I19|. 

In a close vicinity of the critical density ( n c = 0.85 x 10 11 cm -2 for sample Si-2), p{T) becomes non-monotonic, with 
a slow decay after passing through a maximum. The maximum and the subsequent decay are explained as transition 
to the non-degenerate regime and subsequent increasing the scattering rate due to increase in the available phase 
space for the scattered electrons [20]. The non-degenerate regime is however beyond the scope of out study. 

Following the ancient maxim “The simple is the seal of the true”, we conclude on the primary role of the scat¬ 
tering rate (i.e. two channel scattering) rather than scattering time (i.e., two conductivity channels) in this “high- 
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FIG. 3: (Color online) Zero field temperature dependencies of (a) p(T), and (b) u{T). Sample Si-2. 


temperature” transport in the two-phase metallic-like state at zero field. 

In accord with above arguments, we focus further on the p(T) dependence (shown in Fig. 1 of the main text and also 
in Fig. 3 here). On the generic p(T) experimental curve there are two remarkable points: the maxima and inflection 
points, both behave critically and vanish to zero at n = n c E3- The temperature of the maximum is much higher 
than Tkink and is close to the nondegeneracy, whereas the inflection point Ti n fl -C Tp appears to be very close to Tkink- 
The inflection and large p(T) variation is known to be inherent to high mobility samples [2ll423j : as mobility decreases 
(and, hence, critical density increases, and interaction weakens) the magnitude of the p(T) variation diminishes. We 
again conclude that this feature (inflection) is related with electron-electron interactions. 

To summarize, all three notable temperatures, TkinkW) Ti n fl(n), and T dM / dn (i) are close to each other, (ii) behave 
critically, and (iii) are intrinsic to clean samples/strong interactions regime solely. These arguments support our 
suggestion that they indicate a novel “high energy” scale T*(n ) in the two-phase strongly correlated state. 


Comparison of the zero field transport with theory 


The interaction corrections theory, being the low-energy theory, does not include an additional energy scale T*. 
On the other hand, within the framework of the temperature dependent screening model, the large p(T) raise, its 
further saturation, and decay (for low densities) are explained by smearing of the Fermi distribution function in the 
scattering time averaged in the ener gy s pace, that produces negative dp/dT , i.e. by a transition from the degenerate 
to nondegenerate regime above Tp [20j . Within the screening theory the low temperature analytical expansion of 
resistivity is given by [rfllill : 



P 



f T \ 

+ C3/2 (l>) ~ B * 



where all prefactors are positive, 


<?! = 


1 + ^?’ 
Qof 


C-3/2 ~ 


2.646 


Bn = 


7 r 2 p(p + 1) 

6 ’ 


(5) 


Qo = qTF/Zkp oc n -1 / 2 , and / = /(2/cp) is the quasi-2D subband form-factor at the wave vector 2 Rf (/ = 1) in with 
the strictly 2D limit. 

The inflection in this analytical expansion arises due to the third term. In the T/Tp <C 1 limit, the inflection point 
T*/Tp ~ (1 — 0.12-\/n/10 11 ) [25.126 1. however decreases , rather than grows with density. From more exact numerical 
calculations, in Fig. 3b of Ref. 27J, the inflection point also decreases from T*/Tp = 0.7 to T/Tp = 0.2 as density 
grows from 1 to 10 x 10 n cm -2 . Therefore, the calculated T*(n) dependence does not agree with the experimental data 
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shown in Fig. 3 of the main text, where the inflection point increases with density from T /Tp = 0.09 to T/Tp = 0.7 
in the same range of densities. In the latest version of the screening theory [27j, all theoretical p(T) curves have a 


positive curvature with no inflection point (ignoring the weak localization corrections which are beyond the scope of 
our paper and the above theories). 

We believe therefore, that currently no theory capture the main feature of all experimental data for Si-MOSFETs, 
in the “high temperature” correlated regime, including the large p(T) variations and inflection in p(T) whose tem¬ 
perature grows with density. In the absence of a thorough microscopic theory, and having three clear prompts on the 
experimental side, we are eligible to suggest an empirical phenomenological form, that appears to fit the data with a 
minimal set of parameters and links the transport features in a unified picture. The model attaches importance to 
the existence of a novel energy scale T*, and presumes the existence of two parallel scattering channels. Surely, the 
microscopic ground of the novel energy scale requires separate study, that is beyond the framework of our paper. 


Testing the empirical model 

Figure 4 of the main text shows comparison of the calculated p(T) and a a with experimental data. The empirical 
dual resistivity model fits sclf-consistently both, the zero field resistivity p(T) and the magnetoconductivity. The 
parameters of the fit are summarized in the Table below. The factor ft is an order of magnitude smaller than £, 
therefore, the second term in Eq. (3) of the main text becomes important only at high temperatures. The slope, a, 
is almost constant, confirming the linear A (n) dependence (Eq. (1) of the main text). 


TABLE I: Summary of fitting parameters, corresponding to Fig. 4 and Eq. (3) of the main text, pi and pr> = are in 
(kfi/IU), density is in units of 10 11 cm -2 , n c = 0.88, a is in K/10 n cm~ 2 . 


n 

Pd 

Pi 

a 

P (K/T 2 ) 

£ (K 2 /T 2 ) 

1.5 

1268 

14362 

4.53 

-0.0160 

-0.08 

1.996 

901 

9564 

4.35 

-0.0080 

-0.09 

2.5 

662.2 

6937 

4.28 

-0.0043 

-0.11 

3.25 

501.5 

5202 

4.24 

-0.0019 

-0.15 

5.252 

336.14 

3456.6 

4.18 

-0.0005 

-0.19 


On the role of phonons 

In 3D metals any residual weak temperature dependence in p(T) originates from phonon scattering which produces 
the Bloch-Gruneisen behavior, p{T) = po + piT 5 , where the temperature- independent contribution po arises from 
short-range disorder scattering and the temperature dependence (the second term) - from phonon scattering. By 
contrast, the temperature dependent transport in 2D metallic systems at low temperatures, besides weak-localization 
effects, is dominated mostly by electron-impurity scattering dressed with electron-electron interaction effects (or on 
the complementary language - by screened disorder scattering with temperature dependent screening). 

The interaction effects in transport are proportional to (TV) and in order to diminish them and to highlight the 
effect of phonons, we presented in Fig. 0] the resistivity data for low mobility Si-MOS sample (r is smaller by a 
factor of 10 than for the high mobility samples studied in the paper). One can see that below about 2K, logarithmic 
quantum corrections dominate (both, WL and interaction corrections) (2S| (see Fig. QJi). For higher temperatures, 
up to the Fermi energy (dashed curve), ballistic interaction corrections (or temperature dependent screening) take 
over and cause p(T) growing which flattens and then saturates as T approaches Tp, due to nondegeneracy effects 
[20j . For temperatures higher than 100- 200 K, resistivity again starts growing, now due to electron-acoustic-phonon 
scattering. The monotonic p oc T dependence is a consequence of the amount of phonons excited at a given T. In 
GaAs heterostructures, due to effective piezo-coupling the phonon scattering is rather strong [2(1129|- For Si, the 
phonon scattering contribution to the overall scattering rate is much lower, because of the weaker electron-phonon 
coupling mechanism (that is the deformation potential for Si). 

To conclude, it is well-known that phonon scattering in Si-structures contributes essentially to the transport only 
in the vicinity of room temperature, and is irrelevant to the low-temperature transport. Both, the nondegeneracy and 
phonon scattering are irrelevant to the inflection in p(T) which happens at much lower temperatures than the onset 
of phonon scattering. Nevertheless, to be on a safe side, in the main text we analyze our data (kink in d 2 a/dB 2 ) and 
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Temperature (K) Temperature (K) 


FIG. 4: (Color online) Zero field temperature dependencies of p(T) for a low mobility sample in a wide range of temperatures, 
(a) double-log, and (b) linear scale. Sample Si-46. Density: 1 16.5, 2 - 17.6, 3 - 19.8, 4 - 22, 5 - 24.2, 6 - 26.4, 7 - 30.8, and 

8 - 36.3 xlO n cm -2 . Dashed line depicts Tf for various densities. 


inflection (in p(T)) only in the temperature range (i) well below Ep and below the p(T) maximum in order to avoid 
the nondegeneracy effects, and (ii) always below 20K for the explored densities, where the phonon contribution to the 
resistivity in Si-MOSFETs can be neglected with 1% or better accuracy. 


[1] G. Zala, B. N. Narozhny, I. L. Aleiner, Phys. Rev. B 65, 020201 (2001). 

[2] G. Zala, B. N. Narozhny, I. L. Aleiner, Phys. Rev. B 64, 214204 (2001). 

[3] N.Teneh, A.Yu. Kuntsevich, V. M. Pudalov, and M. Reznikov, Phys. Rev.Lett. 109, 226403 (2012). 

[4] V.M.Pudalov, Pis’ma v ZhETF 66, 168 (1997);[JETP Lett. 66, 175 (1997)]. 

[5] Y. Hanein, U. Meirav, D. Shahar, C. C. Li, D. C. Tsui, Hadas Strikman, Phys. Rev. Lett. 80(6), 1288 (1998). DOI: 
10.1103/PhysRevLett.80.1288. 

[6] S. J. Papadakis and M. Shayegan, Phys. Rev. B 57, R15068 (1998). 

[7] X. P. A. Gao, A.P. Mills, Jr., A.P. Ramirez, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 88, 166803 (2002). 

[8] G. Brunthaler, A. Prinz, G. Pillwein, P.E. Lindelof, J. Ahopelto, Physica E 13(2-4), Pages: 691-694 (2002). Arxiv: 

cond-mat/0207170 

[9] S. Ilani, A. Yacoby, D. Mahalu, H. Shtrikman, Science 292, 1354 (2001) 

[10] B. Spivak, Phys. Rev. B 67, 125205 (2003). B. Spivak, and S. A. Kivelson, Phys. Rev. B 70, 155114 (2004). 

[11] Y.V. Stadnik and O. P. Sushkov, Phys. Rev. B 88, 125402 (2013). 

[12] M.W. C. Dharma-wardana, and Francois Perrot Phys. Rev. Lett. 90, 136601 (2003). 

[13] Junren Shi, X. C. Xie, Phys. Rev. Lett. 88, 086401 (2002). 

[14] B. N. Narozhny, I. L. Aleiner, and A. I. Larkin, Phys. Rev. B 62, 14 898 (2000). 

[15] J. W. Clark, V. A. Khodel, and M. V. Zverev, Phys. Rev. B 71, 012401 (2005). 

[16] B.L. Altshuler, D.L. Maslov, and V.M. Pudalov, Physica E, 9, #2, 209-225 (2001). 

[17] D.A. Knyazev, O.E. Omelyanovskii, V.M. Pudalov, and I.S. Burmistrov, Phys. Rev. Lett. 100, 046405 (2008). 

[18] S.V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, V.M. Pudalov, and M. D’lorio, Phys. Rev. B 50, 8039 (1994). 

[19] S.V.Kravchenko, W.E.Mason, G.E.Bowker, J.E.Furneaux, V.M.Pudalov, M.D’Iorio, Phys. Rev. B 51, 7038 (1995). 

[20] S. Das Sarma and E.H. Hwang, Phys. Rev. B 61, 7838 (2000). 

[21] V.M. Pudalov, G. Brunthaler, A. Prinz, G. Bauer, Metal-insulator transition in two dimension, Physica E 3, 79-88 (1998). 

[22] V.M.Pudalov, Chapter in “The Electron Liquid Paradigm in Condensed Matter Physics”, eds. G.F. Giuliani, and G. 

Vignale (IOS Press Amsterdam, 2004), pp.335-351. ArXiv:cond-mat/0405315 

[23] V. M. Pudalov, G. Brunthaler, A. Prinz, and G. Bauer, Pis’ma ZhETF 68, 415 (1998). [JETP Lett. 68, 442 (1998)]. 

[24] A. Gold, V.T Dolgopolov, Phys. Rev. B 33, 1076 (1986). 

[25] S. Das Sarma and E. H. Hwang, Phys. Rev. B 68, 195315 (2003). 

[26] S. Das Sarma and E. H. Hwang, Phys. Rev. B 69, 195305 (2004). 

[27] S. Das Sarma, E. H. Hwang, Qiuzi Li, Phys. Rev. B 88, 155310 (2013). 

[28] B. L. Altshuler, D. L. Maslov, and V. M. Pudalov, phys. stat. sol. (b) 218, 193 (2000). 

[29] X. Zhou, B. Schmidt, L. W. Engel, G. Gervais, L. N. Pfeiffer, K. W. West, and S. Das Sarma, Phys. Rev. B 85, 041310(R) 

















7 


( 2012 ). 


